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Contrary to the conventional empirical law, band gap of dilute nitride semiconductors decreases with nitrogen 
concentration. In spite of a number of investigations, origin of this “large band gap bowing”, is still under 
debate. In order to elucidate this phenomenon, we investigate change of band edge energies of GaN x Asi_ x due 
to nitrogen by using the perturbation theory. It is found that energy shift of conduction band edge is arising from 
mixing between T- and L-states and/or T- and X-states induced by displacement of Ga atoms around N. We 
also found that the valence band edge state shows upward shift in spite of negative potential of nitrogen. These 
results are well understood from character of the wavefunctions and symmetry of the perturbation potential. 


I. INTRODUCTION 

III-V compound semiconductors containing small amount 
of nitrogen atoms are attracting much attention. JMh Exten¬ 
sive studies have been carried out since these materials are 
promising candidate for novel optical devices such as infrared 
lasers and high efficiency solar cells. MM In addition to 
the prospect for device applications, their properties which 
are largely different from other semiconductors, also evoke 
much interest. lflll - [l3ll In particular, large band gap bowing is 
a problem still under debate. Usually, band gap of a mixed 
compound is well described as a linear interpolation of band 
gaps of the constituent materials. Deviation from this empir¬ 
ical law, which is called band gap bowing, is sometimes ob¬ 
served but not so large. However, contrary to the empirical 
law, band gap of some compounds such as GaNAs,GaInNAs 
and GaNP containing small amount of N becomes smaller 
with the N concentration. There have been various models 
and theories to explain this phenomenon. lff4l - [l6ll Analysis on 
the basis of tight-binding theory lfl7l[l8tl . empirical pseudopo¬ 
tential method li20l - l23ll and first principle calculations ll24[ l25ll 
have been carried out. However, its mechanism has not been 
clarified, yet. 

In this paper, to elucidate the mechanism of the large band 
gap bowing, we carry out analysis on behavior of the band 
edges for a typical dilute nitride semiconductor GaN v Asi_ x 
(GaNAs). For this purpose, we carry out perturbation calcu¬ 
lations regarding effects due to nitrogen as the perturbation 
along with bulk GaAs wavefunctions as bases. 

In the next section, first, we show the model used in this 
study. After evaluating matrix elements of the perturbation 
potential, reduction of band gap is calculated. Then, by mak¬ 
ing qualitative interpretation in terms of symmetry of the per¬ 
turbation potential and character of the wavefunctions, we 
show how the reduction of band gap due to nitrogen occurs. 


II. PERTURBATION ANALYSIS 

A. Model 

We consider an NjxNdXNd supercell of GaAs in which one 
of As atoms is replaced by a N atom. Note that this supercell 


consists of AN/ unit cells of the zinc blende structure, that is, 
nitrogen concentration is given by x = 1 /(AN/). Introduction 
of the N atom gives rise to change in crystalline potential. We 
consider following three factors: (i) shift of atomic potential 
from that of As to N, (ii) displacement of Ga atoms adjacent 
to the N atom, and (iii) displacement of As atoms on the sec¬ 
ond neighboring sites to the N atom. Then, the perturbation 
potential is given by 

AV(r) = AV & (r) + AV (ii \r) + AV^\r), (1) 

with 

AV (i} (r) = [V N (r - R,) - V As (r - R,)] , (2a) 

A V (ii) (r) = £ [ y Oa(r + r - Rj - &) - V G /r + r - Rj )], 

(2b) 

and 

AV im (r) = [v As (r - Ry ~ ~Vas (r - Rf)\ ■ (2c) 

f 

where the superscripts (i)-(iii) correspond to the three factors 
noted above. That is, AV (n (r) denotes the potential shift from 
that of As to N, and AV (ll) (r) and AV (, "\r) stem from dis¬ 
placement of atoms adjacent to the nitrogen. In these equa¬ 
tions, Vn, Kas, and V Ga are atomic potentials of N, As, and 
Ga, respectively. For the atomic potentials, we used empir¬ 
ical pseudopotentials presented in Refs. [H HI The vec¬ 
tor R, (Rj) indicates location of a zinc blende unit cell, and 
r = (a/4, a/4, a/4) (a is the lattice constant) is a vector that 
indicates position of a Ga atom within a unti cell. Displace¬ 
ment of atoms adjacent to the nitrogen are written as and 
rjj’ for the first neighboring Ga and the second neighboring 
As, respectively. The indices j and / run through so that 
R j-T+f j and Ry + ijy indicate the positions of the first neigh¬ 
boring four Ga atoms and the second neighboring twelve As 
atoms, respectively. We set jjj so that the Ga atoms approach 
toward the N atom by 0.38 A. Similarly, ijy was determined 
so that the second neighboring As atoms approach toward the 
N by 0.1 A. These values were obtained from total energies 
evaluated by the first principle calculations using the CASTEP 
package. 12^| 
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FIG. 1. Energy band of GaAs calculated by using empirical pseu¬ 
dopotentials. Inset: Reciprocal lattice vectors to be connected with 
one another by the perturbation potential of the present model for 
4x4x4 supercell. Note that equivalent points are excluded in the 
following calculations though they are plotted in the figures. 


B. Perturbation matrices and band gap reduction 

In FigureQ] we show energy band of GaAs calculated using 
the empirical pseudopotentials along with plane wave basis 
functions. Using the wavefunctions of bulk GaAs obtained so 
far, we calculate a matrix element defined by 

dV n k,n'k' = (3) 

where if/„k{r) is a wavefunction of bulk GaAs with a wavevec- 
tor k. The index n denotes the band as shown in Fig. Q] where 
the highest valence band and the lowest conduction band are 
labeled n — 3,4 and n = 5, respectively. 

In the inset of Fig. |T] we show k-states that are connected 
with each other by the perturbation potential plotted on the 
first Brillouin zone of the zinc blende structure for the case of 
N d = 4. These k-states discretely distribute as 

27r 

k = —— (n x ,n y ,n z ), (4) 

N d a 

with n x , n y , and n- integers between -N d and N d . This can 
be proven as follows: Since the perturbation potential has 
translational symmetry with a period N d a in all the x-, y-, 
and ^-directions, the matrix element {fp n k\AV\4i n 'k') must be 
fixed when AV(r) is replaced by AV(r + R) with a lattice vec¬ 
tor of the supercell R. From this and the Bloch’s theorem, 
the wavevectors k and k' in eq. Q must satisfy a relation 
k — kf — (2 n)/(N d a)(n x ,n y ,n z ). In the present case, since 
k' is F, k is given by eq. 0. Note that some points on 
the border are equivalent. For example, (2n/a)(l, 0,0) and 
(27r/a)(-l, 0,0) are identical. Hence, one of them must be ex¬ 
cluded to avoid double counting, though both are plotted in 
the figure. Excluding such equivalent points, we have 4N d 3 
k-points in the first Brillouin zone to be mixed due to the per¬ 
turbation potential; there are 256 points for N d — 4 and 2048 


points for N d = 8 required for the calculations. These k-points 
are the points that are folded onto the F-point of the Brillouin 
zone of the N d x N d x N d supercell. 
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FIG. 2. Upper panel: Absolute values of perturbation matrix ele¬ 
ments between the conduction band edge state (5T) and band states 
are plotted as functions of wave vector k along the L-T-X line. Lower 
panel: Square of matrix elements divided by energy separation are 
plotted along the L-T-X line. The dots on the x-axis denote k-states 
relevant to the calculations for N d = 4. 
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FIG. 3. Components of the matrix element. The curve labeled (i) 
—(iii) show matrix elements for the factor (i)-(fii) given by eqs. (2a)- 
(2c). 

In the upper panel of Fig. [2] we show |Vsr^Jtl. he., matrix 
elements between the conduction band edge state (F-state of 
the 5th band) and various states. Note that matrix elements are 
basicall negative quantities reflecting deep nitrogen potential, 
though we plot absolute values since they are complex quan¬ 
tities. The solid- and dotted- curves show the values for the 
bands n = 5 and n = 6, respectively, plotted as functions of 
k along the L-F-X symmetric line of the zinc blende Brillouin 
zone. The dots on the horizontal axis show the k-states that 
are relevant to the perturbation calculations, as we have shown 
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FIG. 4. Upper panel: Absolute values of perturbation matrix ele¬ 
ments between the conduction band edge state (5r) and band states 
are plotted as functions of wave vector k along the L-r-X line. Lower 
panel: Square of matrix elements over energy separation are plotted 
along the L-f-X line. 


in Fig. [3 

In the lower panel, we plot the quantity IVsr.niP/lesr - 
by the solid- and dotted- curves for the band n — 5 and n = 6, 
respectively. Within the viewpoint of the perturbation theory, 
magnitude of this quantity corresponds to energy change due 
to the perturbation. These curves show that contribution from 
states around the X- and L-points of the lowest conduction 
band is much larger than that from other states. As shown 
in the upper panel, the matrix elements for the 5th and 6th 
bands have large values in a wide area of the Brillouin zone. 
However, except for the X- and L-states of the 5th band, large 
energy separation reduces this quantity. 

To examine further the character of the matrix elements for 
the lowest conduction band n = 5, we show contributions 
from the factors (i)-(iii) described in the previous section in 
Fig. [3] In this figure, the curves labeled (i)-(iii) show compo¬ 
nents of the matrix elements due to the factors (i)-(iii), respec¬ 
tively. Note that, since absolute values of complex quantities 
are plotted, a summation of these components does not make 
the total value. Around the T-point, the factor (i) the N po¬ 
tential takes the largest value of the three. Whereas near the 
X-point, the factor (ii) displacement of Ga atoms is the largest. 
At the L-point, the factors (i) and (ii) are comparable though 
the latter is slightly larger. Contribution from the factor (iii) 
displacement of second neighboring As atoms is smaller than 
others in all the region. As we have shown in Fig. [2] contribu¬ 
tions from the X- and L-states are much larger than that from 
other state. This means that the effect due to the displacement 
of adjacent Ga atom is larger than that due to nitrogen poten¬ 
tial. We will discuss this point later in terms of symmetry of 
wavefunctions. 

In Fig. 0 we plot the perturbation matrix elements for the 
valence states. Similar to Fig. 0 IW.nJtl and |V 4 r>il 2 /|£ 4 r-£;iA:l 
are plotted in the upper and lower panels, respectively. We 
note that the degeneracy of the valence states causes a com¬ 


plex situation because any linear combinations between the 
degenerated states satisfy the Schrodinger equation. This 
gives rise to ambiguity of numerical values. To avoid such 
a situation, we made linear combinations between the degen¬ 
erated states so that the states become proper bases of the ir¬ 
reducible representation of the group they belong. In other 
words, we set the states to have proper symmetry. In the cal¬ 
culations, we set 1 // 4 Y to have x-symmetry. Similarly, we set 
ij/ 4 k, i/j\k and i// 2 k to have x-, y-, and z-symmetry, respectively, 
in drawing Fig. [4] 

Different from the conduction states, contribution from 
states in vicinity of the F-point is large. Thus, the quantity 
|V4r,4i| 2 /|£4r - £ 4*1 could be very large when k is close to the 
F-point. This is because of small denominator arising from 
heavy mass of the hole states. However, the states in vicinity 
of the F are irrelevant to the energy change unless N concen¬ 
tration is very low. As indicated on the horizontal axis for the 
case of Nd = 4, the states that are relevant to energy change 
distribute discretely. Therefore, the states close to the F do not 
affect the energy shift of the band edge state. When nitrogen 
concentration is very low, we have to take these states into ac¬ 
count because when Nd is large the ^-states denoted by the 
dots become denser. However, in this case, as we denote be¬ 
low, the factor 1 /(4A,/ 3 ) added to take nitrogen concentration 
into account reduces the contribution from the states close to 
the T. 
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FIG. 5. Filled circles show calculated band gap as a function of ni¬ 
trogen concentration.. For comparison, theoretical data from Ref. [6] 
are also plotted by crosses. Inset shows energy shifts of the conduc¬ 
tion and valence band edge states. 
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Using the perturbation matrix elements, we evaluate energy 
shift of band edges. For the conduction band edge, the sim¬ 
ple perturbation expansion series converges slowly. This is 
because, as shown in Fig. [2] contribution from the X- and In¬ 
states is larger than that from the T-state. This means that 
higher order terms is larger than the first order term of the 
perturbation series. Therefore, we evaluated the band edge 
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energy by diagonalizing the Hamiltonian matrix given by 


^ ' ' k,:l' , 


k*nk tin.n' 3k,k' "t" 


dVnk.n'k' 

4 N d 3 


(5) 


We fixed as n = n' = 5 because interband effect is negligible 
as we have shown in Fig. [2] The factor 1 /(4/V,/’) is necessary 
to treat N concentration properly. As we have noted, since the 
Nj x N,/ x N,/ supercell contains 4/V,/ : ' of zinc blende unit cells, 
the N concentration is given by x = 1/(4A^ 3 ). 

For the valence band edge, we evaluated energy shift using 
the secnond order formula 


AE V = 


4 N d 


r V4r,4r + 


( 4 N/y 


,3y2 Zj 


W\T,nk\ 


z 


«=? ^ 4 


^4r “ £«£ 


( 6 ) 


As we have seen in Fig. [TJ perturbation matrix elements for 
the valence states are not so large as those for the conduction 
states. As a result, the second order formula yields reasonable 
results close to the value obtained by diagonalizing perturba¬ 
tion matrix. 

In Fig. 0 we plot calculated band gap by filled circles as a 
function of nitrogen concentration. The crosses denote results 
of first principle calculation extracted from Fig. 6 of Ref.fl]. 
Reasonable agreement between results of present theory and 
the first principle calculation indicates that the perturbation 
analysis is valid for dilute nitrides as long as N concentration 
is smaller than 1 %. In the inset, we show energy shift of 
the conduction band edge and valence band edge by filled and 
open circles, respectively. In addition to the lowering of the 
conduction band edge, upward shift of the valence band edge 
state is observed. This seems peculiar since the perturbation 
potential due to nitrogen is negative. 


C. Character of wavefunctions and matrix elements 

In Fig. [ 6 ] we show wavefunctions of GaAs at some symme¬ 
try points. From the top to bottom, wavefunctions at the F-, 
X-, and L- points of the lowest conduction band and the wave- 
function of the valence band edge are plotted along the ( 111 ) 
direction. The As atom locates at the position 0, and a Ga 
atom without displacement locates at +0.25 as indicated by 
open and shaded circle, respectively. For each curve, we also 
attached schematic illustration to show character of the wave- 
function, where local wavefunctions of As and Ga are plotted 
by solid- and dotted- curves, respectively. The wavefunction 
of the r state in the conduction band, t//^\ (r), has peaks at 
the atom positions, indicating ,v- 1 i ke character around an atom. 
This wavefunction consists of anti-bonding coupling between 
an .s , -like orbital of As and an ,v-l i ke orbital of Ga. We also see 
that wavefunction around As is localized whereas that around 
Ga is somewhat extended. In the similar way, we see that the 
wavefunction of the X-state, 1 ^ 5 x(r), consists of anti-bonding 
coupling between an ,v-1 i ke orbital of As and a p-like orbital 
of Ga that has a node at the Ga position. The L-state is some¬ 
what complicated. We see that the wavefunction 1 /r 5 i (r) has 
peaks at the positions deviated from the atom positions. This 
means that, as schematically shown in the inset, the s- and 
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FIG. 6. Wavefunctions of GaAs. the T, X, and L points are plotted 
along the (111) direction. Schematic illustrations to show the char¬ 
acter of the wavefunction are also shown in the insets. The arrows 
show positions of As (N) and Ga. The horizontal dotted lines show 
zero for each wavefunction. 


p-orbitals are mixed within an atom. The large amplitude in 
the middle of the atoms indicates that this wavefunction has 
bonding character. The wavefunction has nodes at the 

atom position, meaning p-like character around an atom. We 
observe that this wavefunction consists of bonding coupling 
between p-like orbitals of As and Ga. 

In Fig. [7] we also plot the perturbation potential along the 
(111) line. For each panel, we added schematic atomic poten¬ 
tials with and without nitrogen by dotted- and solid-curves, 
respectively, to illustrate schematically how the perturbation 
potential arises. Since the N potential is deeper than As po¬ 
tential, the perturbation potential at the As (N) site is negative. 
Due to the spherical character of the pseudopotential, the per¬ 
turbation potential has approximately spherical ,v -1 i ke symme¬ 
try. On the other hand, the perturbation potential around Ga, 
that arises from displacement of Ga atoms, is more compli¬ 
cated. As illustrated in the inset, the displacement of Ga gives 
rise to potential increase on one side of the Ga position and 
decrease on another side. Therefore, the perturbation poten¬ 
tial is nearly anti-symmetric around the Ga position, that is, 
/r-like symmetry along the ( 111 ) direction, although the po¬ 
tential curve seems complicated due to complexity of pseu¬ 
dopotential that is not a simple function. 

Comparison between the wavefunctions shown in Fig. [5] 
and the perturbation potential in Fig. [ 6 ] enables us to make 
interpretation for the perturbation matrix elements shown in 
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FIG. 7. Perturbation potential along the (111) direction is plotted 
by solid curve. The arrows show positions of As (N) and Ga. The 
horizontal dotted line shows zero. Schematic illustration of poten¬ 
tials with and without nitrogen is also drawn to show character of the 
perturbation. 

Fig. 130 

First, we see the matrix element at the T state of the con¬ 
duction band, /fVsr.sr As shown in Fig. [3] contribution of the 
Nitrogen potential is dominant to this first order term /IVsr.sr- 
This is because the wavefunction has large amplitude at the 
As (N) position. The fact that displacement of Ga has small 
effect to -dV 5 r, 5 r is well understood; As we have noted, po¬ 
tential around Ga has p-like symmetry and wavefunction is of 
s-like symmetry. Therefore integrating It/'srl 1 2 3 4 x AV (ll) , which 
has an odd-like character around Ga, results in small value of 
the matrix element. 

Different from the I -state, the matrix element at the X-state 
effect from Ga is the largest as we have shown in Fig. [3] This 
is also understood that the wavefunction of the X-state around 
the Ga site is of p-like symmetry. For the coupling between 
the T- and X-states AVswsx, shift of Ga atoms is essential. 
i/< 5 r has s-like symmetry around the Ga atom, whereas both 
i// 5 x and ,4V /In l have p-like symmetry as seen in Fig. [4] There¬ 
fore, we foresee that multiplication of these three quantities 
becomes an even function around Ga, which enlarges the ma¬ 
trix element AVsr, 5 x- 

As we have noted, the wavefunction 1 J /51 has character in¬ 
termediate between the F- and X- states; i/O/, has both ,y- and 
p-characters. This is also well corresponding to the result that 


AVst, 5 l arises from both the factor (i) and (ii). 

We noted that contribution from position shift of the sec¬ 
ond neighboring As atoms is small. This is also explained 
from symmetry. Although displacement of As atoms is about 
1/3 of that of Ga atoms, contribution might be large because 
of larger number (i.e., 12) of neighboring As atoms. As we 
have noted, atom’s position change leads to perturbation po¬ 
tential with p-like symmetry. As seen from Fig. [4] both the 
F-state and X-state have s-like charge distribution around As 
atoms. From a simple consideration on symmetry, we see that 
the perturbation matrix element for shift of As with k and k' 
F or X has a small value. These results indicate that mixing 
between F and X or between F and L induced by lattice dis¬ 
tortion around N is the main factor of the band gap reduction 
of GaNAs, rather than the nitrogen potential. 

As for the valence state, as shown in Fig.0 the valence edge 
state shows upward shift. This is innegligible, since the energy 
change is as large as 1/3 of that of the conduction state, though 
the valence state is seldom referred. It seems strange that the 
valence state shifts upward since the perturbation potential is 
strongly negative at the N position. As shown in the figure, 
the valence state has p-like symmetry, its amplitude is small 
at the position where the perturbation potential is negative. 
The wavefunction has large amplitude at the peripheral where 
the perturbation is positive. This makes the first order term 
positive. In addition, the second order term of eq. ([6| pushes 
up the valence top state. 

III. CONCLUSION 

In order to clarify the mechanism of the large band gap 
bowing of dilute nitride semiconductors, we carried out anal¬ 
ysis based on the perturbation theory. It was found that reduc¬ 
tion of the conduction band occurs due to intervalley mixing 
between the F state and the X- and L- states. The interval¬ 
ley mixing is induced by displacement of Ga atoms around 
the nitrogen, rather than nitrogen potential itself. This is rea¬ 
sonably understood considering symmetry of wavefunctions. 
The displacement of Ga causes potential with p-like symme¬ 
try, which causes strong mixing between the F state and X- 
(L-) state, since the former has s-like symmetry and the latter 
has p-like symmetry at the Ga position. We also found inneg¬ 
ligible upward shift of the valence top states. Although the 
upward shift in the negative perturbation potential of nitrogen 
seems peculiar, this is also understood from the character of 
the wavefunction. 
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